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We present a detailed analysis of the dynamics of photon transport in waveguiding systems in 
the presence of a two-level system. In these systems, quantum interference effects generate a strong 
effective optical nonlinearity on the few-photon level. We clarify the relevant physical mechanisms 
through an appropriate quantum many-body approach. Based on this, we demonstrate that a single- 
particle photon-atom bound state with an energy outside the band can be excited via multi-particle 
scattering processes. We further show that these trapping effects are robust and, therefore, will be 
useful for the control of photon entanglement in solid-state based quantum optical systems. 

PACS numbers: 42.50.Nn, 42.50.-p, 42.50.Gy, 42.70.Qs 



Over the past years, the conception and development 
of solid-state based quantum optical functional elements 
have received steadily increasing interest [ll-0|. As com- 
pared to other approaches, solid-state-based systems of- 
fer an obvious scalability and handling advantage of the 
resulting devices as well as the utilization of modified 
light-matter interactions through judicious designs of the 
corresponding waveguides' dispersion relations and/or 
mode profiles. 

However, since high-quality samples such as coupled- 
optical-resonator- waveguide arrays (CROWs) (3, [H| have 
become available only recently, there is limited theoret- 
ical work regarding the potential of utilizing modified 
light-matter interaction in (effectively) low-dimensional 
quantum-optical systems. The basic underlying problem, 
i.e., that of a system with discrete levels that is coupled 
to a continuum of states has attracted attention for a 
long time For single photons, quantum interference 
effects in one-dimensional waveguides with an embed- 
ded quantum impurity allow the realization of effective 
energy-dependent mirrors For two or more pho- 

tons, this system induces an effective photon-photon in- 
teraction and even bound photon-photon states that may 
be exploited for efficient control of photon-entanglement 



la-Uj. Except for our work on the one-photon case 



9[, all of the above calculations have been carried out 
in the stationary regime. In particular, the more chal- 
lenging few-photon case has been addressed with sophis- 
ticated Bethe-Ansatz 12, [l^ and Lehmann-Symanzik- 
Zimmermann reduction techniques [l3 | that allow one to 
determine the corresponding scattering matrices for such 
systems. However, these field-theoretical approaches em- 
ploy linearized dispersion relations without band edges. 



In the present Letter, we apply our computational 
framework of time-domain simulations using Krylov- 
subspace-based operator-exponential methods |To| to 
the case of few-photon transport through a quantum im- 
purity in a one-dimensional waveguiding system similar 
to wave packet dynamics in electronic systems (lol . 11 



This allows us to analyze the scattering of two or more 
photons at the quantum impurity in a very general way. 
In particular, for a cosine-type dispersion relation, we are 
able to confirm the existence of two bound photon-atom 
states IJ] . Furthermore, we show how these states can be 
excited and controlled through the photon-nonlinearity 
that is induced by the quantum impurity. This eluci- 
dates the mechanism through which the quantum impu- 
rity can be utilized for controlling photon-entanglement. 
In the field-theoretical approaches discussed above 
14l |. the photon-atom bound states are (due to the ab- 
sence of band edges) energetically shifted to infinity and 
are thus removed from the physically accessible Hilbert 
space. 



Starting from the well-known Dickc-Hamiltonian jl5| . 
we can derive a tight-binding Hamiltonian that describes 
photon propagation in an effectively one-dimensional 
waveguide with cosine-type dispersion relation (such as 
the CROWs of Refs. HQ) that is coupled to a quantum 
impurity as [§] 



x=l 



(1) 



Here, a|, and Ux denote, respectively, bosonic (photon) 
creation and annihilation operators at lattice site x and J 



2 



denotes the corresponding hopping element. The quan- 
tum impurity is modeled as a two-level system (TLS) 
with transition frequency = ^/h that is located at 
lattice site Xg and couples with a coupling element V 
to the modes of the photonic band. When measuring 
energies from the center of the band, the correspond- 
ing dispersion relation is tvjjk = — 2Jcos(fca), where a 
denotes the lattice constant and k stands for a wave 
number that lies within the first Brillouin zone. Finally, 
the TLS is described through the Pauli-operators and 
= o'a; ± iay. 

While being physically intuitive, the above Hamilto- 
nian ([l} does not allow for the most transparent discus- 
sion of the underlying physics. Instead, we find it most 
useful to reformulate the problem in terms of the Hamil- 
tonian 

N-l 

H = (^ala^+i + al^-^^a^^ + Qb^ b 

x=l 

+V {a,,M + alj) + Ub^b {b^b - 1) , (2) 

where we have replaced the TLS by an additional bosonic 
lattice site. More precisely, we have replaced the Pauli- 
operators of the TLS by appropriate combinations of 
bosonic creation and annihilation operators, b^ and 6. 
The ground and excited states of the TLS correspond, 
respectively, to none and a single boson on this addi- 
tional site (TLS site). Unphysical multiple occupancies 
of the TLS site have been addressed through the addition 
of the last term on the r.h.s. of ([2]). This term ensures 
that once the TLS site is occupied, i.e., the TLS is in its 
excited state, adding a further boson to the TLS site re- 
quires the energy U > 0. Thus, Hamiltonians ([T]) and ^ 
are equivalent in the limit U oo and this is the only 
case we consider in this work. The [/-term induces in- 
elastic scattering that allows us to discuss the physically 
relevant processes. For actual numerical calculations, we 
use Hamiltonian ([T]). 

With this reformulation several issues become appar- 
ent. Quantum interference processes associated with the 
coupling between TLS and the waveguide modes induce 
an effective interaction between photons as described by 
the nonlinear term UWb {b'^b — l). While this effective 
few-photon optical nonlinearity is spatially localized to 
the immediate vicinity of the TLS site, this system nev- 
ertheless represents a true quantum-mechanical many- 
particle problem. For instance, Hamiltonian ([2]) looks 
very similar to a bosonic version of the celebrated single- 
impurity Anderson model [l6l | that describes magnetic 
impurities in metals. Therefore, it is suggestive to apply 
methods that have been developed for correlated quan- 
tum systems to the Hamiltonians ([T|) and ([2]) 0, [l2| - 
From Hamiltonian ([2]), it becomes apparent that 
the TLS will induce correlations between two or more 
photons. This raises the fascinating question to what 
extent the TLS can be utilized to engineer this entangle- 



ment and what role the photon-atom bound states play 
in this (note that photon-atom bound states have been 
discussed in a different context before [l7l|). 

To address this question, we have to go beyond station- 
ary calculations that determine the scattering matrices 
of photons in plane wave states for linearized dispersion 
relations where the photon-atom bound states are physi- 
cally inaccessible jl2l4l4|. To do so, we employ our com- 
putational framework which we have described in detail 
elsewhere Q. This framework allows us to analyze both 
the dynamics of multi-photon wave packets that interact 
with the TLS and the dynamics of the TLS itself. Fur- 
thermore, it takes into account all aspects introduced by 
the finite-bandwidth dispersion relation. First, we would 
like to note that on energetic grounds a single photon 
cannot excite the photon-atom bound states described 
above and, therefore, these states are of no relevance in 
single- photon scattering calculations from a TLS d, Q . 
In other words, the TLS (partially) absorbs an incoming 
single photon and a decomposition of the system's initial 
state into the (polaritonic) single-particle eigenstates of 
the Hamiltonian ([T]) does not involve the bound photon- 
atom states. Thus, the excited TLS will eventually de- 
cay into its ground state. However, our reformulated 
Hamiltonian ([2]) suggests that, by virtue of the nonlin- 
ear interaction term, the bound states can, in principle, 
be energetically reached via multi-photon processes. In 
Fig. [TJ we demonstrate that this is indeed possible: A 
two-photon wave packet interacts with the TLS and a siz- 
able fraction of the photon population becomes trapped 
at the TLS site. In other words, once the TLS is ap- 
preciably excited by one of the incoming photons, the 
remaining photon sees a modified (saturated) TLS and 
is thus (partially) scattered into the hitherto unreach- 
able bound photon-atom states via multi-particle scat- 
tering processes. After the scattering is complete, the 
bound photon-atom states are again decoupled from the 
continuum (such as is the case for the scattering of a 
single photon discussed above) and, thus, cannot decay. 
These bound states are of a polaritonic nature, i.e., they 
are multi-moded dressed eigenstates of ^ and ([2]) with 
complex wavenumbers solely induced by the existence of 
the waveguide's finite bandwidth. This implies that a 
fraction of the radiation remains trapped at the TLS site 
in form of a partial occupation of the TLS. 

In order to verify the role of the multi-particle pro- 
cesses, we display in Fig. [2] the time evolution of the 
TLS' excited-state occupation for the scattering of multi- 
photon wave packets with different particle numbers. 
The increase in the trapped photon population with the 
number of photons implies a corresponding increase in 
the rate at which radiation is scattered into the bound 
states. The strength of this interaction further depends 
on the detuning of the TLS relative to the photon fre- 
quency as well as on the strength of the coupling matrix 
element V between TLS and the waveguide modes. In 
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FIG. 1: (Color online) Time evolution (in units of h/ J) of 
transmission (T), reflection (i?), and impurity occupation 
{rih) = (fe^fc) ~ + 1/2) for a two-photon wave packet that 
scatters at a TLS. The TLS has a transition energy 
and couples with coupling strength V = J to the central lat- 
tice site So = 100a of a tight-binding lattice with total extent 
L = 199a and hopping element J. The photons are described 
via boson-symmetric wave packets that are constructed from 
single-particle Gaussian wave functions of width s = 6a with 
wave number k = 37r/4a and initial center Xc = 70a (see 
text and Ref. 0] for further details). All calculations are 
stopped at times not exceeding the transit time, i.e., the time 
the wave packet needs to pass through the waveguide, thus 
avoiding artificial refiections from the system's boundaries. 



Fig. [31 we depict the corresponding dependence of the 
trapped photon population at the TLS for a fixed pho- 
ton wave number k ~ Sir/Aa. Consistent with our above 
interpretation, trapping is most pronounced for zero de- 
tuning 5 ^ Q — hjjjk (recall that huj^^sTr/ia = V^J). Fur- 
thermore, proximity of the TLS resonance frequency to 
the band edge (or cut-off frequency of the waveguide) is 
clearly advantageous for realizing efficient trapping: For 
frequencies near a band edge the multi-particle scatter- 
ing mechanism has to provide less additional energy for 
exciting the energetically closest bound state. Less in- 
tuitive is the fact that there exists an optimal coupling 
strength V^pt ~ J between TLS and waveguide modes 
for which maximal trapping occurs. Wc have confirmed 
these findings for a number of different dispersion rela- 
tions. For instance, we have extended Hamiltonian ([2]) 
to include a next-ncarest-ncighbor hopping term j'^' ^ 
that allows us to significantly modify the cos-type disper- 
sion relation of the tight-binding model (not shown). In 
addition, we have found analogous behavior for strictly 
linear dispersion relations with cut-off at finite energies 
(not shown). 

The above results suggest a certain robustness of the 
trapping effect which we have further analyzed by quali- 
tatively considering losses. This is accomplished by cou- 
pling the TLS to a second waveguide that can de-excite 
the TLS into modes other than those of the original 



FIG. 2: (Color online) Time evolution (in units of ft/ J) of the 
impurity occupation (nt) for initial multi-photon states with 
different photon numbers C that are constructed analogous to 
the two-photon states in Fig. [T] The corresponding system 
parameters are L = 99a, xo — 50a, — and V = J. 

The photon parameters are Xc = 25a, s — 5a, and k = 37r/4a. 
The results for photon numbers C = 3 and C = 4 have been 
obtained with a time-dependent density matrix renormaliza- 
tion group (DMRG) technique as described in Ref. p^ . 




FIG. 3: (Color online) Impurity occupation {nt) in the long- 
time limit (see Fig. [1} after scattering of two-photon states 
for different system parameters V (in units of J) and Q (in 
units of J). The fixed parameters are: L = 199a, xq — 100a, 
Xc ~ 70a, s — 12a, and k = 37r/4a. 



waveguide. The incorporation of this "loss channel" into 
the Hamiltonian ([l} thus proceeds by adding two addi- 
tional terms analogous to, respectively, the first (hopping 
term J') and third term (coupling term V') of the r.h.s. 
of ll]). Clearly, the hopping term J' has to be chosen 
such that the energy of the bound photon-atom states 
that are associated with the first waveguide and the TLS 
alone lies in the band of the second waveguide. In Fig. |4j 
we display the time evolution of of the TLS' excited-state 
occupation for different coupling strengths V' of the TLS 
to such a " broad-band loss waveguide" . The trapping ef- 
fect persists even for rather strong coupling to the loss 
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FIG. 4: (Color online) Time evolution (in units of h/J) of the 
impurity occupation {rib) for initial two-photon states (see 
Fig. [1} where a broad-band loss waveguide has been intro- 
duced. The corresponding system parameters are (loss waveg- 
uide's parameters are primed): L = L' = 399a, xo = x'q = 
200a, n = V2J, V = J, and J' = 2J. The strength of the 
coupling V' (in units of J) to the loss channel is varied. The 
photon parameters are Xc = 175a, s = 6a, and k = Sn/4a. 



FIG. 5: (Color online) Time evolution (in units of h/J) of 
the impurity occupation (nt) for a system where two single- 
photon Gaussian wave packets of width s — 7a with different 
intial positions {x'^^ and aj^^') are launched from different 
sides towards the TLS. The corresponding system parameters 
are L = 199a, xq = 100a, Q = V2J, and V = J. The photon 
parameters are xi^^ = 50a, fc*-^' — 37r/4a = — A;'^-*, and the 
initial position aji^' — 150a -I- A^, is varied. 



channel. If we, for instance, interpret the coupling to the 
loss waveguide as a (admittedly crude) model for fabri- 
cation tolerances that in a quasi-one-dimensional system 
couple strictly guided modes to a continuum of radiative 
modes, we are led to speculate that the trapping effect 
would be observable in experimentally accessible systems. 

Finally, we have analyzed the possibility of tuning the 
trapped photon population at the TLS site. To do so, we 
have prepared two identical single-photon wave packets 
on different sides of the TLS and have launched them 
towards the TLS. By changing their initial relative sep- 
aration to the TLS, we can exert some control over the 
multi-particle processes in the Hamiltonian ([2]) (see Fig. 
[S]). While the dynamcis is (expcctcdly) rather distinct 
for the different cases, we observe monotonic behavior 
of the trapped population: maximal trapping occurs for 
symmetrically launched pulses with zero relative initial 
distance. For increased distances the trapped population 
decreases to zero once there is no overlap of the pulses at 
the TLS site. 

In conclusion, we have analyzed the dynamics of pho- 
ton transport in wavcguiding systems in the presence of a 
TLS within the context of a quantum many-body frame- 
work. Our reformulation ([2]) allows us to identify strong 
multi-particle processes that may be utilized to excite 
and control photon-atom bound states. In turn, this fa- 
cilitates trapping of radiation at the TLS. In addition, 
we have shown that this trapping effect exhibits a cer- 
tain degree of robustness and can be found in a number 
of systems. Since few-photon (or low intensity) coherent 
states are superpositions of a few Fock states only (those 
that we have discussed in the present work), we expect 



that the excitation and control of the photon-atom bound 
states and associated effects will also occur in such situ- 
ations. 

Finally, we would like to emphasize the generality of 
our approach which is capable of treating systems with 
arbitrary dispersion relations and atom-field couphng 
strengths both in real and momentum space. Thus, the 
trapping of the photon population and its control sug- 
gest that such systems may be exploited for engineer- 
ing photon entanglement as well as for the realization of 
quantum logic circuits in a number of systems that range 
from silicon integrated optical elements all the way to su- 
perconducting quantum circuits for microwave photons. 
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